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Abstract 



We consider the renormalisation of the Fayet-Iliopoulos D-term in a 
softly-broken supersymmetric gauge theory with a non-simple gauge group 
containing an abelian factor, and present the associated /5-function through 
three loops. We also include in an appendix the result for several abelian 
factors. We specialise to the case of the minimal supersymmetric standard 
model (MSSM), and investigate the behaviour of the Fayet-Iliopoulos cou- 
pling £ for various boundary conditions at the unification scale. We focus 
particularly on the case of non-standard soft supersymmetry breaking cou- 
plings, for which £ evolves significantly between the unification scale and 
the weak scale. 



I. INTRODUCTION 



In abelian gauge theories with N = 1 supersymmetry there exists a possible invariant 
that is not allowed in the non-abelian case: the Fayet-Iliopoulos D-term, 

L = £ J V(x, 6, 9) d 4 9 = ^D(x). (1.1) 

In previous papers |[| @ we have discussed the renormalisation of £ in the presence of 
the standard soft supersymmetry-breaking terms 

Lsb = -(m 2 )^Vi - Q/^' fc M-0fc + \h ij ^i<t>j + \m\X + h.c.) (1.2) 
The result for f3^ is as follows: 

& = ^£ + & (1-3) 

where 0% is determined by V-tadpole (or in components D-tadpole) graphs, and is inde- 
pendent of £. Although in Refs. QTJ , p| we restricted ourselves to the abelian case, it is 
evident that a D-term can occur with a direct product gauge group (Gi <S> G^ • • •) if there 
is an abelian factor: as is the case for the MSSM. In the MSSM context one may treat 
£ as a free parameter at the weak scale ||, in which case there is no need to know /?£. 
However, if we know £ at gauge unification, for example, then we need f3^ to predict £ at 
low energies. Our purpose in this paper is first of all to give the result for (3% through three 
loops for a general direct product gauge group. For simplicity of exposition, we restrict 
ourselves in the main body of the paper to the case of one abelian factor, postponing 
the more general result (which is complicated by the possibility of "kinetic mixing" Q 
between different abelian factors) to an Appendix. We shall then specialise to the case of 
the MSSM, and perform some running analyses to determine the size of £(Mz) for various 
choices of boundary conditions at the unification scale Mx- 



II. THE GENERAL CASE 



First of all, for completeness and to establish notation, let us recapitulate the standard 
results for the supersymmetric theory. We take an N = 1 supersymmetric gauge theory 
with gauge group H a G a and with superpotential 

6 2 



W($) = -Y^^^k + -//"<I\<Iv (2.1) 
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We will be assuming here that the gauge group has one abelian factor, which we shall 
take to be G\. We shall denote the hypercharge matrix for G\ by y . At one loop we have 



16tt 2 /?W = glQ a = g z a [T(R a ) - 3C(G a )] , (2.2a) 
16tt 2 7 W. = P\ = l -Y^Y okl -2'£g 1 a [C(R a )] i j , (2.2b) 

where R a is the group representation for G a acting on the chiral fields, C(R a ) the cor- 
responding quadratic Casimir and T(R a ) = (r Q ,)~ 1 Tr[C(i? Q ,)] , r a being the dimension 
of G a . For the adjoint representation, C(R a ) = C(G a )I a , where I a is the r a x r a unit 
matrix. Note that T{R X ) = Tr[y 2 }, [C(#i)]% = (^ 2 )V At two loops we have 

(16tt 2 ) 2 ^) = 2g 5 a C(G a )Q a - 2<£r- 1 Tr [PC(R a )} (2.3a) 

(i67r 2 )V 2 ^ = \-Y 3mn Y mvi -2j2glC(R a r j S i n ]P n p + 2 y £gtC(R a ) i j Qa. (2.3b) 



For completeness and later reference, we also quote here the general result for /?^ RED< - 3 \ 
which is a straightforward generalisation of the result of Ref . [H] : 

^dredo) = Sr-iglY^Y^P^CiRaYi + Gr-'gl £ ^Tr [PC{R a )C{Rp)} 

P 

+ r-'glTr [P 2 C(R a )\ - Qr^gl £ Q p gpt [C{R a )C{Rp)\ 

P 

- 4r- l g 5 a C(G a )Tr [PC(R a )) + g 7 a Q a C(G a ) [AC(G a ) - Q a ] , (2.4) 

We recall that gauge anomaly cancellation requires 

Tr[yC(R a )] = (2.5) 

and naturalness (or cancellation of U\ -gravitational anomalies) requires 

Tr[y] = 0. (2.6) 

The diagrams contributing to f3% through three loops for a general non-simple gauge 
group are essentially the same as those depicted for the pure abelian case in Ref. [§], but 
reinterpreting internal gauge and gaugino propagators as ranging over all gauge groups in 
the direct product. Potential new 3-loop graphs (involving a 3-point gauge vertex, or a 
gauge/gaugino vertex) give contributions which vanish due to anomaly cancellation (such 
as C(G a )Tr[yC(R a )]). It is then relatively easy to generalise the abelian result to the 
general case. We find 
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IQn 2 ^ = 2g{Ii \y m 



-4g>iTr ym 7 



2 -v(l) 



(2.7) 
(2.8) 



(167T 



(3)DRED' 
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4Tr [yWP] - -Tr [ytftf 1 



= -6(167r 2 ) 2 Tr [Vt (2) 
+ 2Tr [yP 2 m 2 ] - 24((3) ]T (? 2 Tr 
+ 12C(3) £ (7 2 Tr [3^M:ffC(^) + c.c] 
- 96C(3) £ ^M^T^C^C^)] 



a,j9 



- 24C(3) { ^^M a M^TV[^C(^)C(^)] + c.c. 

a,/3 



(2.9) 



where [[ 



^F 2 m 2 + \m 2 Y 2 + /r%- + 2Y ipq Y i 



jpr 



m 



2\r 



•.-S^jJJIf^C^ (2.10) 



^ = / i ^ + 4$:^M /3 [C( J R /3 )] i 



(2.11) 



with {Y 2 Yj = Y M Y jk i, {h 2 ) l j = h lkl h jk i. These results are computed using the DRED' 
scheme, which is a variant of DRED defined so as to ensure that /3-functions for physical 
couplings have no dependence on the e-scalar mass 0. Most of the terms in Eq. ( |2.9|) 
correspond in a simple way to the analogous terms in Eq. (5.2) of Ref. @; the only subtle 
point being the MM*g 4 terms, where one sees easily that only in the case of Fig. 15(e) 
of can the two gaugino masses belong to different gauge groups (G a ). Thus the last 
term in Eq. ( |2.9| ), and the MM*g A terms from the terms involving H, come entirely from 
this particular figure. 

It was proved in Ref. (Jj in the pure abelian case that if the m 2 dependence in (3^ takes 
the form Tr[m 2 A], then we have 



Tr[yA] = 2%. 

9i 



(2.12) 



It is easy to see that the proof extends to the direct product case, and indeed we can 
check Eq. ( g3g ) explicitly using Eqs. and ( ggp , JEgaD and 
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III. THE MSSM 



We now specialise to the case of the MSSM. The relevant part of the MSSM 
tential is: 

W = H 2 t c Y t Q + H x b c Y h Q + H x r c Y T L 

where Y t , Y b , Y T are 3x3 Yukawa flavour matrices. 
The gauge /^-functions are given at one loop by 

16vrX = b a g 3 a , 



where 



33 
5 ' 



1, 



-3, 



and our U\ coupling normalisation corresponds to the usual one such that g\ 
For the anomalous dimensions of the chiral superfields we have at one loop: 



where 



16^7? 


= P tc = 


167r 2 7£ } 


= P^c = 




= Pq = 


16^7$ 


= P T c = 




= Pl = 


16vr 2 7 ff 


= P Hl = 




— Ph 2 - 








Cfr = 




c Q = 




C T c = 




c L = 



TY t Y} - 
2^ 



YiY b + Y}Y t -2C Q , 
2Y T Y^ — 2C r c j 
Y}Y T - 2C L , 
= Tr [F T tr T + 3YjY b 

2C L 



2C L , 



3Tr 



U + hi 
U + hi 

4 fl 2 , 3 2 , 1 2 
3 </3 ^ 4^2 ^ 60^1' 

4^2 ^ 20^1 - 



At two loops || the anomalous dimensions are given by: 



(i6tt 2 ; 


o (2) 


= -2Y t (P Q + P H2 ) y/ - 2P t cC t c + 2(±b l9 f + |6 35 r|), 


(3.6a) 


(16tt 2 ; 




= -2y (P Q + PffJ F 6 f - 2PcC fe c + 2(i^ x 4 + §m 3 4 ) ; 


(3.6b) 


(i6tt 2 ; 




= -Yj (P t c + P^ 2 ) Y t - Y} (P fec + P Hl ) Y b 








- 2P Q C Q + 2{±hg A x + \b 2 g A 2 + |6 3 ^|), 


(3.6c) 


/ Hi 

(167T 2 ; 


2 7 f c ) 


= -2Y T (P L + PffJ y T t - 2P T cC T c + lb l9 i 


(3.6d) 


(i6^ 2 ; 


2 7? 


= -Y}[P T . + P Hl ]Y T - 2P L C L + ±b ig t + lb 2 gl, 


(3.6e) 


(i6^ 2 ; 




= -3Tr[Y b P Q Y; + Y^P bc Y b ] - Tr[Y T P Q Y} + Y}P T cY T ] 








-2C L P Hl + f- b ig t + lb 2 gl, 


(3.6f) 


(167T 2 ; 


2 tS 


= -3Tr[Y t P Q Yj + YfPrYt] - 2C L P H2 + ^b x g{ + \b 2 g\. 


(3.6g) 



We now turn to the soft couplings. The quantities W and H defined in Eqs. ( |2.1U| ), ( |2.11| ) 
are given by 



-iMM 



iMM 



W t c = (2m 2 tc + ^i\ 2 )Y t Y} + AY t m 2 Q Y t j + 2Y t Y}m\ c + Ah t h\ - 8C$ 
W b , = (2ml + *™ 2 Hl )Y b Yl + AY b m 2 Q Y b ] + 2Y b Y b ] ml + Ah b h\ - 8C>. 
W Q = (m 2 Q + 2m 2 H2 )Y t % + (m 2 Q + 2m 2 Hl )Y*Y b + {Y?Y t + Y h ] Y b ]m 2 Q + 2Y}m 2 t< Y t 

+ 2Y b ] m 2 bc Y b + 2h\h t + 2h\h b - 8C^ M , 
W T c = (2m 2 c + Am 2 Hi )Y T Y} + AY T m 2 L Y} + 2Y T Y}m 2 rC + 4h T h\ - 8C. 



MM 



w L 



[m, 



+ 2m 2 )Y}Y T + 2Y}m 2 T cY T + Y}Y T m\ + 2h\h T - 8C 



MM 
L > 



Tr 



Q(m 2 Hl + m 2 Q )Y b ] Y b + Qm 2 bc Y b Y^ + 2(m 2 Hi + m 2 L )Y}Y T + 2Y}m 2 rC Y T 
+ 6hlh b + 2h\h T ] - 8Cf M , 
W H2 = 6Tr[(m| 2 + m 2 Q )Y t ] Y t + m 2 c Y t Y? + h\h t ] - 8C" L 



iMM 
'L > 



(3.7) 



where 



U t c — 


11^3 


2 , 3 2 + A| Ml | 2 , 2 , 


/-iMM 

(-/ 6 c — 


11^3 




(jMM _ 


11^3 




qMM _ 






qMM _ 


||M 2 


2 9i + l \Mi\ 2 9l 



(3- 



and 



4htYl + 4C£, 



H bc = 4h b Yj + ACj¥, 
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H Q = 2{Y}ht + Y^h b )+4C^, 
H T c = Ah T Y} + 4C£, 
H L = 2Y}hr + 4C?, 
H Hl =Ti[6Y^h b + Y^h T ]+AC M 



L > 



i/^ = 6Tr[y//i t ] + 4Cf, (3.9) 



where 



C£ = lM 3 gi + ±M l9 i, 

C% = iM 3 gl + ±M l9 l 

C% = lM 3 gl + \M 2 g 2 2 + ±M l9 l 

<% = lM l9 l 

Cf = \M 2 g\ + l Q M ig f. (3.10) 

With all these subsidiary definitions we can now give the results for (3^ up to three loops. 
We have 

2/5(1) _ o /3„ TVl™2 2 r. 2 , 2 , ^2 2 , 2 



16tt 2 /3^ ; = 2y / f 5l TV[m| - m£ - 2m 2 c + m\ c + m 2 c - m 2 Hi + m^J (3.11) 



(16tt 2 ) 2 $ 2) = -4 v /|( 7l Tr(m 2 3 P g - m|P L - 2m%P tc + m\ c P h , 

+ m 2 c P rC - m 2 Hl P Hl + m 2 H2 P H2 ) (3.12) 



(IBtt 2 ) 3 ^ = V / !^i|-6(167r 2 ) 2 ^ - 4/?g' - §4" + 2$ 



2\2fl(3) .0(3) 5 /o(3) , ofl( 3 ) 



+ C(3)(-24/4 3) + 12/?g> - 96/3® - 48/3®) I , (3.13) 



where 

/f 

/f 

$ 
41 
4° 
41 



= TV(mJ 7 « - mhf - 2m?„ 7 <? + mfog> + m?, 7 <? - m?,^ + m?, l7 g>), 
= TV(W<,P(, - W L P L - 2H 7 .P t . + HvPj. + W T .P T , - W ffl P ffl + W H2 P„ 2 ), 
= Tv(H 1 Q H Q - H\H L - 2H\,H r + Hl.H t . + Hl.H T . - H^H„ 2 + H^H„,), 
= •n(m%P% - m\Pl - 2m%P 2 + m \,Pl + m 2 ,P?. - m 2 m P 2 m + m 2 H2 P 2 2 ), 
= Tl(W Q C Q - W L C L - 2WrC e + W„.C V + W T .C T < - W m C L + W„ 2 C L ), 
= ^{HqCq" - H L Cl" - 2H,.C^' + H ¥ Cg' + H T ,C™' - H H fi" " + H„ 2 C£") + c.c, 

_ o/piMM/i siMM^y or<MM s~i i n^M r< 1 nMMfi \ 

— 0\yjQ <_/Q — Ui - /Ufc U(C ■+- Ofcc Ofcc -+- O r c W C J 

= 3(|Cj| 2 - |Cf | 2 - 2|Cf | 2 + |Cf | 2 + |C r ^| 2 ). (3.14) 



6 



We shall now present our MSSM results specialised to the commonly considered case 
where only the 3rd generation Yukawa couplings are significant. We also take the gaugino 
masses to be real. Writing At, X b and A T for the 3rd generation couplings, Eq. 
becomes 





— 0\ 2 _ OCir 
^^t 


P, 

r b c 


— 9 x 2 on. 


P^ 


— \ 2 4- \ 2 — OC^ 


P „ 


- 9\ 2 — Or r 


r> 
Ph 


= X T - 2C L 


Pu c 


— —2Cf<: 




= — 2Cb c 


Pr 


= -2C Q 


P e c 


= — 2C T c 


Pn 


= -2C L 


Phi 


= X 2 T + 3X 2 b - 2C, 


Ph 2 


= 3A 2 - 2C L , 



(3.15) 

where {t, b, Q, r, L} now refers to the 3rd generation, and {u, d, R, e, N} refers to either 
of the 1st or 2nd generations. Eq. ( ft. 6a ) now takes the form 



(167r 2 ) 2 # = -2A 2 (P Q + P H2 ) - 2P t .C tc + 2{±b x gl + f^), (3.16a) 
(167r 2 ) 2 7 J 2) = -2A 2 (P Q + P Hl ) - 2P b oC bc + 2(U l9 t + |6 3 ^)> (3-16b) 



(167r 2 ) 2 7 g ) = -A 2 (P ta + P H2 ) - A 2 (P bc + P Hl ) 

- 2P Q C Q + 2(±b l9 t + \b 2 g\ + lb 3 gl), (3.16c) 

(167r 2 ) 2 T ^ } = -2A 2 (P L + P Hl ) - 2P TC C T c + lb l9 t, (3.16d) 

(167r 2 ) 2 7 i 2) = -A 2 [F rC + P Hl ] - 2P L C L + ±b l9 t + lb 2 gt (3.16e) 

(levr 2 ) 2 ^? = -2P uC C tc + 2{±b l9 { + §& 3 <7 3 4 ), (3.16f) 

(levr 2 ) 2 ^ = -2PM. + 2{±b 1 gi + psgt), (3.16g) 

(levr 2 ) 2 ^ = -2P R C Q + 2{±b l9 \ + \b 2 g\ + psgl), (3.16h) 

(16vr 2 ) 2 7 ^ = -2P eC C T c + \bxgi (3.161) 

(IGtt 2 ) 2 ^ = -2P N C L + p l9 t + \b 2 g\ (3.16j) 

^2\2 (2) _ o\2r c , , r> 1 \2r 



{^Yfui = -*K[Pq + *V] - K\Pl + PA 

- 2C L P Hl + p l9 * + \b 2 g\, (3.16k) 
(167r 2 ) 2 7 g = -3A 2 [F Q + P tc ] - 2C L P H2 + U l9 t + lb 2 gl (3.161) 
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Correspondingly, we retain only the three 3rd generation trilinear soft couplings h t 
h t = A b X b and h T = A T X T . Eq. (|3.7| ) now becomes 



4A 
4A 



2/ 2 i 2 , 2 
77V + f^Q + 171 



h 2 + A 2 t ) 



bK m b , 



+ m 2 Q +m 2 Hl +A 



4A 



,777 



2 TC + m 2 Hl +ml + A 2 T )-8C. 



2X 2 T (m 2 Hi + m 2 T c + m\ + A 2 ) 



W Q = 2X 2 (m 2 c +m 2 Q + m 2 H2 +A 

w L 

W R 

w N 



MM 



-8Q 
— 8C MM 
+ 2A^( 



ml? + + n^Hi + A 



gQMM 



MM 



8Cf M 



— 8C MM 
—8C MM 
—8C MM 
—8C MM 



]C MM 



w, 



6A 



m 



2 Hi + m\ c + m 2 Q + Al) + 2X 2 T (m 2 Hi + m\ + m 2 rC + A 2 ) 



gQMM 
L 



h 2 = 6X 2 (m 
Eq. ( p.9| ) now becomes 



2,2,2 
H2 + 77V + m Q 



At 



IC MM_ 



H t c 


= 4AA? 


+ 4Cf , 


H b c 


= 4A b X 2 b 


+ <, 


Hq 


= 2(A t X 2 t + A b X 2 b ) + 4Cj 


H T c 


= AA T X 2 T + 4:C^, 


H L 


= 2A T X 2 


+ 4Cf 




= 4C£, 




Hd c 


= <, 




Hr 


= 




H e c 


= 4C r A f, 




H N 


= 4Cf , 




Hhi 


= 6A b X 2 b 


+ 2A T X 2 + 4Cf 


Hh 2 


= 6A t X 2 t 


+ 4Cf . 



Eqs. ( J3.ll ), (3.12) now become 



+ m 2 c 



9i ( m Q + 2t77^ — m\ — 2m 2 N — 2m 2 c — ^m 2 uC + m 2 bc + 2m 2 , c 
+ 2m? c 



777 



77i 



+ 777 



H 2 h 
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(16vr 2 ) 2 A (2) 



c - -4^1 (m^Fg + 2m 2 R P R - m 2 P L - 2m 2 N P N - 2m\ c P t , - Am 2 uC P u c 
+ m 2 bc P b c + 2m 2 dc P d c + m 2 rC P T c + 2m 2 eC P e c - m 2 Hi P Hl + m 2 H2 P H2 ) . (3.20) 



Finally, Eq. ( pj.l4|) is replaced by 



,( 2 ) 



(2) 



2„(2) 



.,(2) 



.,(2) 



Aml^S +ml^ ] +2m 2 dcl fJ 



+ m 2 rC ^) + 2m\ c ^) 



2 (2) . 2 (2) 



4 3) = W Q P Q + 2W R P R - W L P L - 2W N P N - 2W t cP t c - m u oP u o + W b .P b . + 2W d .P d . 
+ W TC P TC + 2W eC P eC - W Hl P Hl + W H2 P H2 



Pi 



(3) 
'6 



rj-2 I o fx 2 t_t2 r> tt2 r> t_t2 j tt2 , tt2 , r> rr2 
rr2 , c) tt2 tt2 , r_r2 



4 



(3) 

«4 



m 2 P 2 + 2mi P 



R 1 R 

+ m 2 TC P 2 c + 2m 2 c P 2 



4m 2 c P 2 c + m 2 c P 2 +2m 2 c P 2 



m 2 H P 2 Hl + m 2 H2 P 2 H2 



4 3) = WqCq + 21W Q - W/ L C L - 2WjvC l - 2WtcC t c - 4^C, + W*C* + 2W d *C b c 
+ W TC C T c + 2W e cC T c - W Hl C L + W H2 C L 

iq( 3 ) _ of u n M i ofj n M tj n M oh n^ 1 oh n M a h n M _i_ h n M _i_ oh n M 
P£ 6 — Z{HqUq + ZiIrUq — Hl^l ~ Zsin^l ~~ ^-n t cU t c — 4il u cU t c + n. b ^ hc + Zn d cU bc 

i tt riM i o H S~iM tt riM , tt /~iM\ 



4 



(3) 
«8 



I (J M 1 2 _ I (7 M 1 2 
I Q I I _L I 



2|C5| 2 + |C£T + |G 



M|2 



yM|2 



(3.21) 



IV. RUNNING ANALYSIS 



As we mentioned in the Introduction, if we have no prejudice as to the value of £ 
at the gauge unification scale Mx, then we may as well treat £ as a free parameter at 
the weak scale 0, and the running of £ is irrelevant. However it is conceivable that the 
underlying theory at scales beyond Mx may favour certain values of £(Mx), and then 
the running of £ would need to be considered. We shall see that for currently popular 
choices of boundary conditions at M^-namely, the minimal supergravity scenario, and 
the Anomaly Mediated Supersymmetry Breakdown (AMSB) scenario-the running of £ is 
determined predominantly by the first term on the RHS of Eq. (|1.3|) between Mx and 
Mz, and hence to a good approximation we have 

i(M z ) « ^4^\aMx). (4.1) 
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For instance, we find from Eqs. (|2.7| ), ( |2.8| ) that universal soft masses at Mx imply 
P^\m x ) = (Mx) = 0, using Eq. ( |2.6|) and the fact that it follows immediately from 
Eq. Q2.2b|) using gauge invariance and anomaly cancellation (Eq. (|2.5| )) that 

Tr[^7 (1) ] = 0. (4.2) 

Moreover, it is easy to show, using the result for (3^1 from Ref. @, that if we work 
consistently at one loop, then Tr[ym 2 ] is scale invariant. So if initially £ = Tr[3^m 2 ] = 0, 
then £ remains zero under (one loop) RG evolution. With typical universal conditions at 
Mx with soft masses m and M ~ m , A ~ m , we find (using three loops for (3^ and two 
loops for the other /3-functions) that £ ~ 0.001m 2 , at Mz- 

Another favoured set of boundary conditions is those derived from anomaly mediated 
symmetry breaking (AMSB) 0. Here the soft masses are given by 



(m 2 )^|m3|V^ ; (4.3) 
where rri3 is the gravitino mass. In fact, since the AMSB result is RG invariant, it applies 

2 

at all scales between M x and M Z - We then find from Eqs. ( |2.7| ), ( |2.8| ) that up to two 
loops, we may write 

16n 2 ^ = gi \ms lV^Tr[^( 7 - T 2 )]- (4-4) 



Gauge invariance and anomaly cancellation combined with Eqs. fl2.2b| ) and ( |2.3b| ) yield 

i 

Tr[^7 (1) ] = Tr[^( 7 (2) - (7 (1) ) 2 )] = 0, (4.5) 

and so (3^ vanishes through two loops. Therefore to a good approximation ^(Mz) will be 
given by Eq. (11), and once again will be negligible at Mz if it is zero at Mx- 



However, if non-universal scalar masses at Mx are contemplated, then the effects of 
might be significant-as was noted in Ref. [0, for instance. Another context where 



might play a role is that of non-standard soft-supersymmetry breaking [IT]. This is 
because with the non-standard terms (for example <p 2 <p* terms) the result that Tr[ym 2 ] 
is one-loop scale invariant is not preserved. It follows that even with universal boundary 
conditions for m 2 and £ = at Mx, £ becomes non-zero at Mz even with one-loop 
running. In the current context of the MSSM with the 3rd generation dominating, the 
additional soft terms are given by 
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£goft = m^ Hl iP H2 + A t \ t H{Q c + A b X b H* 2 Qb c + A t X t H;Lt c + h.c. 



(4.6) 



Now in Ref. [IT] we assumed, in fact, that £ was zero at Mz] here we explore the more 
natural assumption that £ = at the unification scale. We follow Ref. jLlJ in dropping the 
explicit //-term from the superpotential, since it can be subsumed into . With given 
values at Mx for and for the universal parameters A, M and mo, and for a given tan f3, 
we adjust A t = A b = A T = A (at Mx) to obtain an acceptable electroweak vacuum. As in 
Ref. [JIT]] , we have made allowance for radiative corrections by using the tree Higgs minimi- 
sation conditions, but evaluated at the scale Msusy ~ mo- In Fig- 1 we show (for illustra- 
tive values of M, and A) the region of the mo, tan f3 plane where this can be achieved. 



Allowed region with M=200GeV, A=0, m =150GeV,|(M x )=0 




tan beta 



Fig.l: The region of the mo, tan/3 plane corresponding to an acceptable electroweak vac- 
uum, for M = 200GeV, = 150GeV ; A = and £(My) = 0. The shaded region 
corresponds to one or more sparticle or Higgs masses in violation of current experimental 
bounds. 



For comparison, we show in Fig. 2 the corresponding region 
for £(M S tj S y) = 0. We notice that it is qualitatively similar, though slightly larger. 
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Allowed region with M=200GeV, A=0, m =150GeV,§(M .. J=0 




2 4 6 8 10 12 14 16 18 20 22 24 

tan beta 



Fig. 2: The region of the tuq, tan/? plane corresponding to an acceptable electroweak vac- 
uum, for M = 200GeV, = 150GeV, A = and £(Msusy) = 0. The shaded region 
corresponds to one or more sparticle or Higgs masses in violation of current experimental 
bounds. 



Note that this Figure differs slightly from Fig. 1 of Ref. [11 1. This is because we 



have incorporated one-loop corrections to the Higgs mass and because we have taken 
account of the increasingly stringent experimental bounds (in particular increasing 
at Mx to get acceptable chargino masses). For mo = 640GeV and tan/3 = 8, we 
find A = 1.07(1.01)TeV, A t (M SVSY ) « 661(627)GeV, A t (M SU sy) ~ 664(630)GeV, 
A b (M SVSY ) w 491(469)GeV. (The pairs of numbers correspond to £(M X ) = 0, 
^(M SU sy) = respectively.) The spectra obtained for £(Mx) = and for £(M SU sy) = 
are given in Table 2. We see that there are significant differences, especially in the masses 
of H, A and H ± . On the other hand the chargino and neutralino masses are unaffected, 
with a LSP neutralino. 
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639 


649 


t 2 


319 


334 




604 


615 


bo 


776 


780 


Ti 


625 


639 


To 


663 


658 


^ L 


792 


794 




793 


785 


dr 


796 


798 


d,u 


781 


785 


P T 


664 


657 


P E> 


632 


646 


V 


657 


650 


u e 


659 


652 


h 


116 


116 


H 


231 


291 


A 


230 


290 


E ± 


244 


302 


Al 


120 


120 


A2 


201 


201 


Al 


68 


68 


X2 


116 


116 


X3 


167 


167 


X4 


234 


234 


5 


521 


521 



Table 1: Spectra (in GeY) for £(M X ) = and for £(M SU sy) = 0, with M = 200GeV, 
m = 640GeV 7 A = 0, = 150GeV at M x , and with tan/3 = 8. 

Finally, in Table 2 we give the values of £(Msusy) for some typical points in the allowed 
region of Fig. 1. We see indeed that £(M SU sy) is quite sizeable. 
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m (GeV) 


tan 3 




640 


8 


-5.07 x 10 4 


700 


6 


-5.48 x 10 4 


700 


8 


-5.02 x 10 4 


700 


16 


-5.15 x 10 4 


800 


6 


-5.61 x 10 4 


800 


8 


-4.90 x 10 4 


800 


16 


-4.75 x 10 4 



Table 2: Values for £(M SUSY ) with $(M X ) = and with M = 200GeV, A = and 
= 150GeV at M x . 

We have verified that the same results are obtained if we either (1) Perform the RG 
evolution in the £-uneliminated theory and then eliminate £ (via its equation of motion) 
at low energies or (2) Eliminate £ at M x , and evolve to low energies with the (modified) 
^-eliminated /3-functions. For a general discussion of the equivalence of these procedures, 
see Refs. |J, @. 

V. CONCLUSIONS 

In this paper we have extended the results of Ref. |2| for the renormalisation of the 
Fayet-Iliopoulos D-term to the case of a direct product gauge group, and applied the 
result to the MSSM. With standard soft supersymmetry breaking and universal boundary 
conditions at M x , then £ is negligible at low energies if £(M X ) = 0. However with non- 
standard soft breakings (and/or non-universal boundary conditions for the standard ones) 
we find significant effects even for £(Mx) = 0. In the non-standard breaking case, the 
effect is especially marked for the masses of H, A and H^, which decrease significantly 
when £ is taken into account. 
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APPENDIX A: GENERAL RESULT FOR SEVERAL ABELIAN FACTORS 



In this Appendix we give the general results for the case of a direct product group 
with several abelian factors. As we mentioned earlier, the situation is complicated by 
the possibility of "kinetic mixing" between the different abelian factors. We can 
accommodate this possibility by introducing a matrix of couplings for the abelian factors. 
Suppose that the gauge group is (Ui) a nS= a +i ^<*> where the G a , a = a + 1 . . . n are non- 
abelian. The gauge couplings are then g a p, where g a p = g a 5 a f3, a — a + 1 . . .n, with a 
similar form for the gauge ^-functions. The gaugino masses also form a matrix M Q/ g with 
an analogous structure, as do their /^-functions. Suppose the hypercharges of the abelian 
factors for a given representation are y a , a — 1 . . . a. Then we define 

a 

y a = E y^pa, a = 1 ... a (Al) 
p=i 

and a generalised quadratic Casimir 

a n 

C(R) = Y,y a y a + E 9 2 a C(R a ). (A2) 

a=l a=a+l 

The Fayet-Iliopoulos couplings now form a vector £ a , a = 1 ... a, and we have the matrix 
equation 

^ = g-^ 9 i + h- (A3) 



We can now give the explicit general results, starting with the gauge /3-functions and 
anomalous dimension. At one loop, 

16tt 2 /? 9 « = gQ (A4) 



where 



Q af3 = Tr [y a yp\ , a, (3=1... a 

Qa3 = 9lQaS a l3, « = Q + 1 . . . U 



(A5) 



and 



16tt 2 ( 7 «)\- = = \Y ik % kl - 2C{R) 1 



(A6) 



At two loops, 



15 



(16^ 2 ) 2 ( 7 (2) )*, = - [Y ]mn Y m ^ + 2C{Rf 3 5\] P% + 2 (Q^y^ + g 4 a Q a C(R a 



31 



(AT) 



and 



a, (5 — 1 ... a 



(16tt 2 ) 2 (/3 9(2) ) a/3 = -2g ai Tr [Py,y? 
(lQn 2 ) 2 (P 9{2) ) a =2glC{G a )Q a -2glr- 1 Tx[PC{R a )}, a = a + l...n. 

At three loops we have 

(16tt 2 ) 3 (/9 9DRED(3) ) q/3 = g ai {SY ikm Y jkn P n m (y,y^ J t + 6Tr[P^ /3 C( J R)] 

+ TrfP 2 ^] - 6 £ QkaTi^W^a] 

K,A = 1 

-6 jr gtQ K Tr[y,ypC(R K )]}, a, 0=1... a 

K=a+1 

(>dred(3)^ = 3r- 1 ^y^-y jfen p™ m c( J R a )^ + Gr^Tr [PC(R a )C(R) 



+ r a^a Tr 



F 2 C( J R Q )] - Qr- l gl £ Q /tA Tr[C( J R a )^ A ] 

K,A=1 



(A8) 



Qr-'gl £ ^Q K Tr[C( J R a )C( J R K )]-4r- 1 ^C(G' a )Tr[FC( J R Q )] 



K=a+1 



+ glQ a C{G a ) [4C(G Q ) - g a ] , a = a + 1 . . . n. 



(A9) 



For the Fayet-Iliopoulos couplings we have at one loop 



16tt 2 



(i) 



Tr 



^ Q m 2 



and at two loops 



16tt 2 



= -4Tr 



ct = 1 ... a 



M 1 ) 



Finally, 

(16tt 2 ) 3 ^ ; 



) = -6(16yr 2 ) 2 Tr [y a m 2 7 
+ 2Tr 



(2) 



(A10) 



(All) 



4Tr 



y a wp 



- -Tr 

2 



!V a P 2 m 2 J - 24C(3)Tr p a WC{R) 
+ 12C(3)Tr p^ifC^*^) + c.c] - 96((3)Tr[y a C MM \R)C{R)] 
- 24C(3) {Tr\y a C M (R)C M *(R)} + c.c.} , (A12) 
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where 



-Y 2 m 2 + i 
2 2 



W^- = (^ 2 m 2 + im 2 7 2 + /I 2 )*,- + 2Y tpq Y Jpr (m 2 ) r q - 8C MM * (i?)',, 



C M (i?)= £ M a/3 ^+ E 9 2 a M a C(R a ), 

a,/3=l a=a+l 



C MM (R)= J2(MM*) a/3 y a y p + Yl M a M* a glC{R a ). (A13) 

Q,/3=l a=a+l 



17 



REFERENCES 



[1] I. Jack and D.R.T. Jones, Phys. Lett. B473 (2000) 102 

[2] I. Jack, D.R.T. Jones and S. Parsons, Phys. Rev. D62 (2000) 125022 

[3] A. de Gouvea, A. Friedland and H. Murayama, Phys. Rev. D59 (1999) 095008 

[4] B. Holdom, Phys. Lett. B166 (1986) 196; 

F. del Aguila, G.D. Coughlan and M. Quiros, Nucl. Phys. B307 (1988) 633; 
F. del Aguila, Acta Phys. Polon. B25 (1994) 1317; 

F. del Aguila, M. Cvetic and P. Langacker, Phys. Rev. D52 (1995) 37; 
K.S. Babu, C. Kolda and J. March-Russell, Phys. Rev. D54 (1996) 4635; 
K.R. Dienes, C. Kolda and J. March-Russell, Nucl. Phys. B492 (1997) 104; 
K.S. Babu, C. Kolda and J. March-Russell, Phys. Rev. D57 (1998) 6788 

[5] I. Jack, D.R.T. Jones and C.G. North, Nucl. Phys. B486 (1997) 479 

[6] I. Jack and D.R.T. Jones, Phys. Lett. B333 (1994) 372 

[7] I. Jack, D.R.T. Jones, S.P Martin, M.T. Vaughn and Y. Yamada, Phys. Rev. D50 
(1994) R5481 

[8] J.E. Bjdrkman and D.R.T. Jones, Nucl. Phys. B259 (1985) 533 

[9] L. Randall and R. Sundrum, Nucl. Phys. B557 (1999) 79; 

G. F. Giudice, M.A. Luty, H. Murayama and R. Rattazzi, JHEP 9812 (1998) 27; 
A. Pomarol and R. Rattazzi, JHEP 9905 (1999) 013; 

T. Gherghetta, G.F. Giudice and J.D. Wells, Nucl. Phys. B559 (1999) 27; 
M.A. Luty and R. Rattazzi, JHEP 9911 (1999) 001; 

Z. Chacko, M.A. Luty, I. Maksymyk and E. Ponton, JHEP 0004 (2000) 001; 

E. Katz, Y. Shadmi and Y. Shirman, JHEP 9908 (1999) 015; 
I. Jack and D.R.T. Jones, Phys. Lett. B465 (1999) 148; 

J.L. Feng and T. Moroi, Phys. Rev. D61 (2000) 095004; 
G.D. Kribs, Phys. Rev. D62 (2000) 015008; 
S. Su, Nucl. Phys. B573 (2000) 87; 

J.A. Bagger, T. Moroi and E. Poppitz, JHEP 0004 (2000) 009; 
R. Rattazzi et al, Nucl. Phys. B576 (2000) 3; 

F. E. Paige and J. Wells, |hep-ph/0001249| 

[10] T. Falk, Phys. Lett. B456 (1999) 171 

[11] I. Jack and D.R.T. Jones, Phys. Lett. B457 (1999) 101; Phys. Rev. D61 (2000) 095002 



18 



